arXiv:astro-ph/0512093vl 5 Dec 2005 


Mon. Not. R. Astron. Soc. OOP. [TIT^ ('2002') Printed 5 February 2008 (MN style file v2.2) 


On The Detection of Artefacts in Spectro-Astrometry 


E.Brannigan^*, M. Takami^’^, A. Chrysostomou^, J. Bailey^ 

^ Centre for Astrophysics Research, University of Hertfordshire, Hatfield, HERTS ALIO 9AB, UK 
^Subaru Telescope, 650 North A’ohoku Place, Hilo, HI 96720, USA 
^Anglo-Australian Observatory, PO Box 296, Epping, NSW 1710, Australia 


Accepted 2005 ??. Received ????; in original form ???? 


ABSTRACT 

We demonstrate that artificial bipolar structure can be detected using spectro- 
astrometry when the point spread function (PSF) of a point source suffers distortion in 
a relatively wide slit. Spectro-Astrometry is a technique which allows us to probe the 
spatial structure of astronomical sources on milliarcseond (mas) scales making it pos¬ 
sible to detect close binaries and to study the geometry and kinematics of outflowing 
gas on scales much smaller than the seeing or the diffraction limit of the telescope. It 
is demonstrated that distortion of the PSF, caused by tracking errors of the telescope 
or unstable active optics during an exposure can induce artihcial signals which may 
be misinterpreted as a real spectro-astrometric signal. Using simulations we show that 
these may be minimised by using a narrow slit relative to the seeing. Spectra should 
be obtained at anti-parallel slit position angles (e.g., 0° and 180°) for comparison in 
order to allow artificial signatures to be identified. 

Key words: Spectro-Astrometry Line: Formation Line: Profiles Stars: Pre-Main 
Sequence ISM: Jets and Outflows. 


1 INTRODUCTION 

Spatial resolutions of optical-IR facilities have dramatically 
improved over the last decade. Indeed, the Hubble Space 
Telescope and adaptive optics on ground based 8-10m tele¬ 
scopes have provided spatial resolutions of 0.05 arcseconds, 
leading to a better nnderstanding of the nature of a variety 
of astronomical objects. Despite this, even higher spatial 
resolutions are desired to study the geometry of extrasolar 
planetary systems, populations and formation of close bi¬ 
naries, mechanisms of mass ejection and accretion in young 
and evolved stars and the nature of active galactic nuclei. 
Optical-IR interferometry has begun to provide resolutions 
on milliarcsecond (mas) scales, however its applicability is 
limited due to the complexity of the technique. 

Spectro-astrometry is an alternative approach to study 
milliarcsecond structures at optical-IR wavelengths. The ba¬ 
sic concept is to measure the relative position of the source as 
a function of wavelength using either an intensity weighted 
centroiding algorithm or profile fitting (see Bailey 1998a for 
a more detailed description). The resultant “position spec¬ 
trum” will show structure on any spectral feature that is 
displaced from the centroid of the continuum source: e.g., 
an emission line arising in a binary companion, outflowing 
jets or any other structure that is not perfectly symmetric 
about the centroidal source, or whose displacement from the 
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source varies with wavelength. The smallest spatial scale ob¬ 
served is limited by the seeing or the diffraction limit of the 
telescope, however this technique indeed provides astronom¬ 
ically useful information at spatial scales much smaller than 
these. 

The concept was introduced in the 1980’s as “Differen¬ 
tial Speckle Interferometry” or “Chromatic Position Differ¬ 
ence” (e.g.,Beckers 1982), and demonstrated to detect close 
binaries (e.g., Sorokin & Tokovinin 1985). This broad band 
technique was limited by the necessity to correct for atmo¬ 
spheric dispersion which can be much larger than the ob¬ 
served structure. This difficulty is removed by measuring the 
displacement of a spectral line at high resolution. Aime et 
al. (1988) were able to measure spectro-astrometric struc¬ 
ture using special instrumentation. See Bailey (1998a) for 
review of these works. 

More recently Bailey (1998a) adapted this technique 
using a standard long-slit CCD spectrograph and coined 
the term “spectro-astrometry”. The method is similar to 
that used by Solf & Bohm (1993), who studied a jet from 
a young stellar object (YSO) on subarcsecond/arcsecond 
scales. Bailey (1998a) revised this method and has shown 
that with good spatial uniformity of the CCD coupled with 
a pixel scale that allows excellent sampling of the seing pro¬ 
file positional accuracies as small as Imas are achieved (see 
Takami 2003). Bailey has also shown that the technique is 
a powerful tool for discovering binary companions toward 
young and evolved stars, separating their spectra, and study- 
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ing the kinematics of YSO jets/winds and narrow line re¬ 
gions of active galactic nuclei (AGN). Since this method 
doesn’t require any special instrumentation, its popular¬ 
ity has steadily increased in recent years. Studies made for 
young binaries and YSO jets include Bailey(1998b), Takami 
et.al. (2001,2002,2003), Davis et al. (2001,2003), Whelan, 
Ray & Davis(2004) and Baines et al. (2004). Measurements 
using an integral field unit (IFU) were first conducted by 
Garcia, Thiebaut & Bacon (1999). A similar approach using 
VLTI-AMBER is investigated by Petrov (2003), Marconi et 
al. (2003) and Dougadous et al. (2003). Since the positional 
accuracy achieved depends on the PSF size as well as the 
photon noise, using an interferometer for this approach has 
the potential to achieve information on spatial structures at 
microarcsecond scales (Bailey 1998a). 

Instrumental effects which may compromise spectro- 
astrometric observations have not been fully investigated. 
These effects may indeed create false detections of close bi¬ 
naries or bipolar outflows. Possible instrumental effects in¬ 
clude; misalignment of the spectrum with GCD columns, any 
departure of the CCD pixels from a regular grid, imperfect 
flatfielding or charge transfer deficiencies in the CCD (Bai¬ 
ley 1998b; Takami et al. 2001). To eliminate instrumental 
effects, Bailey (1998a,b) proposed to obtain position spec¬ 
tra at anti-parallel slit position angles (e.g.,0° and 180°) via 
rotation of the instrument. On comparison of these spectra, 
any real signal from the object will change sign, while the 
instrumental effects remain constant. We can thus eliminate 
instrumental effects by subtracting one position spectrum 
from the other. This method is followed by Takami et al. 
(2001,2002,2003) and Baines (2004). Garcia et al. (1999), 
Davis et al. (2001,2003) and Whelan et al. (2004) however 
eliminate instrumental effects by polynomial fitting of the 
continuum position. 

Bailey (1998a) reports a systematic effect upon obser¬ 
vation of sharp unresolved lines. This is caused by either 
telescope tracking errors or unstable active optics during 
the exposure, and appears in the position spectrum as a 
P Cygni type profile. Bailey warns that the target spectral 
lines should be well resolved to observe true positional dis¬ 
placement. This instrumental effect could be eliminated in 
some cases by Bailey’s method, but this cannot be done by 
polynomial fitting of the continuum position. In this paper 
we show that this instrumental effect can appear even if 
the line profile is fully resolved as was the case in position 
spectra of RU Lupi, obtained using the 3.5m New Technol¬ 
ogy Telescope (NTT). Takami et al. had previously detected 
the presence of bipolar outflow in this object using the 3.9m 
Anglo-Australian Telescope (AAT). Indeed the observed po¬ 
sitional displacement in the NTT looks similar to that seen 
in the AAT data. We show that the signal detected in the 
NTT data is not real, motivating us to investigate the effect 
of uneven illumination on spectro-astrometry, how to avoid 
it and whether previously detected bipolar outflows are in 
fact real. 

In §2 we describe the spectro-astrometric observations 
obtained using EMMI on the 3.6m NTT and the results 
in which artificial signatures were detected. These artefacts 
are identified by qualitatively inspecting position spectra ob¬ 
tained using anti-parallel slit position angles. In §3 we (1) 
show that the observed position spectra are explained by 
distortion/motion of the stellar image in a relatively wide 


slit, (2) show that the extent of the artefacts, in a given set 
of observing conditions, can be greatly reduced by using a 
slit width narrow in comparison to the seeing, and (3) show 
an example of the effect these artefacts have on real spectro- 
astrometric data. In §4 we give some recommendations on 
how to monitor the data for the false signatures of bipolar 
structure. In Appendix A, we show that the bipolar outflows 
reported by Takami et al. (2001,2003) are real and not due 
to the instrumental effect described here. 


2 OBSERVATIONS AND RESULTS 

Medium resolution spectroscopy of RU Lup was obtained 
with the 3.5m NTT in the early morning of April 8 
2004, when the seeing was measured to be 0.8”. The ob¬ 
servations were made using the REad Medium Disper¬ 
sion (REMD) mode of the ESO Multi-Mode Instrument 
(EMMI) with the #6 grating and an arcsecond slit. The 
wavelength coverage spanned 6100 to 6800A with resolv¬ 
ing power = 5000. The spectra were obtained at 4 position 
angles: 0°, 90°, 180°, 270° by rotating the instrument. HeAr 
lamp spectra were used to calibrate the wavelength scale and 
a pixel scale of 0.165” with the MIT detector provided good 
sampling of the seeing profile. Flat fields were obtained by 
combining many exposures of the spectrograph illuminated 
by a halogen lamp. 

The data were reduced using the FIGARO package pro¬ 
vided by Starlink, following Bailey (1998a,b) and Takami 
(2001,2003). After subtracting the bias level and flat- 
fielding, “position spectra” were determined by fitting the 
seeing profile at each wavelength with a gaussian function. 
Global instrumental effects in the position spectra are re¬ 
moved by polynomial fitting. This allows us to see the posi¬ 
tional displacement of any emission/absorption features at 
each slit angle relative to the continuum. In addition to the 
position spectra, the intensity spectra were obtained by sub¬ 
tracting the bias, flat-fielding, subtracting the adjacent sky 
and extracting bright columns on the GCD. 

The spectrum of RU Lup shows a wealth of emis¬ 
sion lines as observed by Lago & Pension (1982) and 
Takami (2001). The intensity and position spectra of the 
Ha A6563A, Hel A6678A & [SII] A6731A emission lines and 
Lil A6708A absorption line, obtained using EMMI-NTT, are 
shown in Figure 1. The Ha emission line shows a large offset 
in the position spectrum, ~25mas from the continuum posi¬ 
tion and the blue and red shifted wings are displaced in dif¬ 
ferent directions. Similar displacement in Ha is reported in 
Takami (2001). Such structure usually indicates the presence 
of a bipolar outflow (see Appendix A). However a qualitative 
inspection of position spectra, obtained using anti-parallel 
slits (i.e. 0°&180°) in this data, reveals that the displace¬ 
ment seen in anti-parallel position spectra have equal sign 
indicating that the signature must be systematic. In addi¬ 
tion, the displacement observed with slit position angle 0° 
is less than that observed with position angle 180°. In fact 
two position spectra observed at a position angle of 0° ex¬ 
hibit dissimilar profiles indicating that the effect responsible 
for the spurious displacement signal is time variable and as 
such cannot be removed even by obtaining position spectra 
at anti-parallel slit angles. 

In addition to Ha A6563A and Hel A6678A, evidence for 
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Figure 1. Intensity and position spectra of RU Lup obtained using EMMI-NTT. From left to right; Top: Ho A6563A, Hel A6678A 
&: [SII] A6731A emission and Lil A6708A absorption lines. The intensities shown are arbitrary; Lower four panels: spectro-astrometric 
signal (mas) for each corresponding spectral feature detected at 4 position angles 0°, 180°, 90°&:270° in the spectrum of RU Lup. The 
bottom two panels show position spectra for anti-parallel slit angles 90° &: 270° in each of the spectral features, similarly in the third 
and fourth panels from the bottom position spectra for anti-parallel slit angles 0° &: 180° in each of the spectral features is shown. In the 
blue wing of [SII] A6731A emission feature the angular displacement is seen to change sign in corresponding anti-parallel slits, whereas 
displacement near zero velocity of the [SII] A6731A line (marked with a dot-dashed line), in the permitted emission lines and the Lil 
A6708A absorption line does not change sign in anti-parallel slits. 



Figure 2. Image capture in slit viewer 7/4/2004 using 3.5m NTT. 
Pixel scale: 0.4". Left and Centre: distorted image of the Classical 
T Tauri Star (CTTS) SZ68. Right: distorted image of the CTTS 
SZ69 


bipolar structure is seen in other bright permitted lines in- 
eluding Fell A6238/6248/6456/6516A. The Lil A6708A ab¬ 
sorption line exhibits displacement similar to that of Ha and 
the bright permitted emission lines, although in any given 
spectrum the direction of displacement in Li I absorption is 
of opposite sign to that of the emission lines. For example, 
at a position angle of 180° the blue and redshifted compo¬ 
nents of the bright permitted emission lines show positive 
and negative displacement respectively, whereas the Li I ab¬ 
sorption line shows negative and positive displacement in 
the blue and red wings respectively. 

Conversely, the direction of displacement due to 


blueshifted forbidden line emission reverses sign at anti¬ 
parallel slit positions (0° vs. 180° ; 90° vs. 270°). This is best 
seen in the [SII] A6731A emission line in Figure 1. However 
a large, sharp displacement is observed close to zero veloc¬ 
ity in the position spectrum at both [SII] A6731A and [01]^ 
A6300A emission lines. This displacement is observed in all 
the position spectra for [SII] A6731A (marked in Figure 1 
with a dot-dashed line) and [01] A6300A. This sharp dis¬ 
placement does not reverse sign in anti-parallel slit position 
angles and therefore this is likely to be due to instrumental 
effects and not to actual signal from the object. 


3 INDUCED SIGNAL DUE TO A DISTORTED 
STELLAR IMAGE IN A WIDE SLIT 

The instrumental effects described in §2 were not observed 
in previous observing runs at the AAT described in Takami 
et al. (2001,2003). There are a number of key differences be¬ 
tween these earlier observations and those reported in §2. 
Primarily, the typical seeing conditions differ dramatically 
(0.8" and 1.5"/2.0") while a wide slit of 1.0" was used in 

^ not shown as the effect is in this line is the same as that in 
the [SII] line, but the data at the wavelength of [OI] is of lower 
quality 
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Figure 3. Top left: Intensity spectrum of RU Lup obtained on 8/4/2004 using EMMI on the 3.6m NTT. Lower left: simulated position 
spectra for a distorted PSF using 3 slit widths, 0.85", 0.5" Sz 0.17", flux ratio Fi : Fi = 0.1 and two point sources located off centre of the 
slit with identical spectra. Lower right: geometry used in the simulations to show that a distorted PSF can reproduce well the artefacts 
seen in the data. Top right: Close up view of the Hel, [SII] and Lil emission and absorption lines.Lower top right: positional displacement 
due to Hel, [SII] and Lil obtained from the simulations, where Fi : Fi =0.1 and slit width = 0.5". The simulated displacement is similar 
to that seen in the data. 


all of the observations. These artefacts are seen in observa¬ 
tions when the seeing is smaller than the slit width used. 
Secondly, the stellar image suffered intermittent distortion 
during the later observations, as shown in Figure 2, due to 
an unstable active optics system. In §3.1 it is explained how 
these combined effects can induce an artificial signal in the 
position spectrum, and the extent of the artefacts that can 
be expected to occur for a given set of observing conditions is 
investigated. In §3.2 an example of the effect these artefacts 
can have on real spectro-astrometric signal is simulated. 

3.1 Concepts and Simulations 

When the intensity distribution from a given source is non- 
uniform across the slit the resultant spectrum can be slightly 
blue or redshifted. This blue/redshifted emission due to “un¬ 
even illumination” or “inaccurate centering”, (e.g., Bacciotti 
et al. 2002; Ardila et al. 2002) arises from an offset of the 
light source from the center of the slit, resulting in a devi¬ 
ation of the angle of incidence according to da = dy/fcol, 
where fcoi is the focal length of the collimator. According to 
the grating equation: 

mA = (sin a ± sin/?) (1) 

where m is the diffracted order, A the diffracted wavelength, 
a is the angle of incidence wrt. the normal and (5 the an¬ 
gle of diffraction wrt. the normal, this leads to a subsequent 


deviation in the angle of diffraction which in turn is trans¬ 
lated to a shift in the spectrum. Therefore, the shift in the 
spectrum, dA, is proportional to the position of the light 
source from the center of the slit axis dy, dA = const.x 
dy. This constant value can be readily determined from the 
spectral resolution of the instrument. For the REMD mode 
of EMMI-NTT a 1.0" slit provides R= 5000, corresponding 
to AA = I. 2 A giving us the relation: 

dA [A] = 1.2 X dy [arcsec] (2) 

As such, if a stellar image is distorted or elongated in the 
slit, perhaps due to tracking errors of the telescope or un¬ 
stable active optics, the spectrum close to the edges of the 
slit is slightly blue or redshifted. This geometry would be 
replicated by a bipolar outflow located within the slit. In 
order to determine the affect of these artefacts on spectro- 
astrometric observations simplified simulations of two point 
sources were carried out with their positions as depicted in 
Figure 3. Given that the point sources have spectra Fi/(A) 
and F 2 /(A), where Fi and F 2 are constant and /(A) is a 
normalised spectrum, the spectra are recorded at the detec¬ 
tor as Fi/(A -I- AA) and F 2 /(A — AA). Thus the centroidal 
position at each wavelength is given by: 

y ,,, Ax Fi/(A + 1.2Ay) - F 2 /(A - 1.2Ay) 

^ 2 Fi/(A-bl.2Ay)-bF2/(A-1.2A3/) 
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Parameter 

Description 

Distortion 

Binary 

Binary Sz Distortion 

PS 

Pixel Scale 

0.17 

0.17 

0.17 

^min 

start wavelength 

6100 

6100 

6100 

^max 

end wavelength 

6850 

6850 

6850 

PAbin 

Binary position angle 

- 

45° 

45° 

ASbin 

Binary separation 

- 

0.1" 

0.1" 

Fi 

Flux of primary binary component 

1.0 

3.0 

3.0 

F2 

Flux of secondary binary component 

- 

1.0 

1.0 

^distort 

Percentage of flux in distortion 

0.1 - 0.9 

- 

0.1 - 0.9 

^^distort 

Distortion position angle 

135° 

- 

135° 

ASdistort 

Distance between distortion and “real image” 

0.1" - 1.0" 

- 

0.1" - 1.0" 

FWHM 

Gaussian simulated seeing 

0.2" - 1.7" 

0.2" - 1.7" 

0.2" - 1.7" 

Aw 

Slit width 

0.1" - 1.4" 

0.1" - 1.4" 

0.1" - 1.4" 


Table 1. Parameters used in the simulations. 


To confirm that the artefacts seen in this data are a 
result of uneven illumination of the slit simulations of a dis¬ 
torted PSF were carried out and position spectra obtained 
using equation and the geometry described in Figure 
3. The simulations reproduce a number of the artefacts ob¬ 
served in the data (Figure 1), specifically (see Figure 3): (1) 
evidence for bipolar structure exists in the bright permit¬ 
ted emission lines; (2) displacement seen at the Li I A6708A 
absorption line has opposite sign to that of the emission 
lines in a given spectrum; and (3) sharp displacement close 
to zero velocity is exhibited in the forbidden emission lines 
[01] A6300A and [SII] A6731A. The extent of the displace¬ 
ment varies with the ratio Fi:F 2 , and in fact the relative 
displacement in the red and blue wings in the emission lines 
in a given position spectrum also depends upon this ratio. 
The extent of the displacement dramatically decreases as 
the slit width becomes smaller (see Figure 3). 

The point source geometry depicted in Figure 3 does 
not provide a good representation of observing conditions, 
as such the seeing conditions are simulated by convolving 
the two point sources by use of a gaussian function. The 
simulated intensity distribution is then described as: 

F{x,y) = Gi{x,y)+ G 2 {x,y) (4) 


where: 

Gi{x,y) 

G2{x,y) 


2.773 

{FWHMY-k 

2.773 

{FWHMYn 


— 2 773 + 

Pj g (FWHM)'2 


— 2 773 + 

P 2 e (FWHM)'-‘ 


where FWHM is the full width half maximum of the seeing 
profile, (A xi, A yi) is the position of the primary point 
source, (A X 2 , A y 2 ) the position of the secondary point 
source, Fi and F 2 are the flux of primary and secondary 
point sources respectively. Combining equations 0 & 0 
results in: 


, , _ J Jx F{x,y) /(A -f AX{y)) da: dy 
// F(a;,j/)/(A-|-AA(t/))da;d?/ 


The two point sources, having being allocated the rel¬ 
evant relative flux, are convolved by our simplified model 
of the seeing and subsequently shifted according to equa¬ 
tion 1^ A position spectrum is then obtained as described 
by equation 13 

A simulated position spectrum (Figure 4), obtained as 
described above, shows displacement on a scale smaller than 
those artefacts obtained by the unconvolved point source. 



Wavelength (A) 


Figure 4. Position spectrum obtained as described in §3.1, where 
the parameters slit width, seeing AS distort and F distort are 1.0'', 
0.8", 0.7" and 0.5 respectively. 


This is because much of the light falls close to and on the 
center of the slit and as such suffers less of a “shift”. In this 
instance a 1.0" slit is simulated with FWHM of 0.8", as was 
the case for those observations shown in Figure 1. The two 
point sources were defined 0.7" apart, with equal flux (i.e. 
Fi:F 2 = 1) and both off center of the slit, as described in 
Figure 3. The angular scale of the displacement simulated 
for Ha emission (~ 25mas) compares well with the artefacts 
observed in the position spectra obtained at the NTT. 

In order to investigate the limitations imposed on 
spectro-astrometry due to systematic effects, the extent of 
the artefacts that can be expected when uneven illumination 
of the slit occurs (perhaps due to tracking errors or instru¬ 
mental concerns) for a given set of conditions is simulated. A 
single point source is allocated at the centre of the slit with 
flux Fi and subsequently convolved by use of a gaussian 
function. A distortion in the resulting PSF is simulated by 
“relocating” a percentage of the flux attributed to the point 
source (Fdistort) to a new position off centre of the slit. The 
results of these simulations are shown in Figure 6, in which 
the full extent of the artificial signal, measured from peak 
to trough, is shown for a range of seeing, slit width, distance 
between the “real image” and distortion (ASdistort) and the 
amount of flux attributed to the distortion (Fdistort), the val¬ 
ues of which are given in Table 1. The results are shown for 
an Fdistort of 0.1 and 1.0 for given slit widths and ASdistort- 
From this we can clearly see that for small ASdistort where 
the distortion in the PSF cannot be detected by eye artifi¬ 
cial spectro-astrometric signal is still generated on a scale 
00 5mas in Ha, above the typical detection limit of Imas 
for previous spectro-astrometric observations (Takami 2001, 
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Seeing (arcsec) 



Figure 5. Simulated angular displacement in Ha, measured from peak to trough, for a range of seeing, slit width, distance between the 
“real image” and distortion {^Sdistort) the amount of flux attributed to the distortion {Fdistorts which gives the ratio Fi:F 2 ), the 
values of which are given in Table 1. Top: Artificial signal generated for a given slit width (solid lines) for a range of seeing values, with 
Fi:F 2=1.0 (left) and Fi:F2=0.1 (right). Bottom: Artificial signal generated for a given separation between “real image” and distortion, 
/^Sdistort (solid lines) for a range of seeing values, with Fi:F2=1.0 (left) and Fi:F2=0.1 (right). 


2003). In fact when the input parameters are similar to the 
conditions during earlier observations described in Takami 
2001, 2003 (i.e. seeing 1.5"coupled with a slit width of 1.0") 
the displacement in Ha shows an angular scale of ~7 - 8mas, 
similarly above the typical detection limit. Therefore it is im¬ 
perative to obtain data at anti-parallel slit angles in order to 
determine, by qualitative inspection, if any detected signal 
is in fact from a real structure in the object, even when the 
PSF appears symmetric. 

The model suggests that the effect can be minimised 
by defocusing the telescope. For example, Figure 5 shows 
that given a slit width of 1.0" an increase of 50% in the 
seeing parameter (from 0.8" to 1.2") results in a reduction 
of the angular scale of the simulated artefacts by a factor of 
2. However, the use of a slit width narrow in comparison to 
the seeing reduces the extent of the artefacts to a greater 
extent. From Figure 5 it can be seen that given a seeing of 


0.8" a decrease in the slit width by 50% (from 1.0" to 0.5") 
results in a reduction of the angular scale of the simulated 
artefacts by a factor of 3. As such it is clear that a reduction 
in the slit width is the best option in order to minimise the 
magnitude of these artefacts. 


3.2 Binary Simulations 

In order to determine (1) what, if any, effect the use of a 
relatively narrow slit might have on the detection of any 
real spectro-astrometric signal from the object, and (2) how 
this systematic effect will affect real signals from the object, 
the convolved model was adapted with the introduction of 
an artificial binary companion. The intensity spectrum al¬ 
located to the secondary binary component is the same as 
that in the primary with line to continuum ratio increased 
by a factor of 2 in order to allow a “real” spectro-astrometric 


Artificial Signal (mas) Artificial Signal (i 
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Figure 6. Binary simulations with (right) and without (left) a distorted PSF. The true displacement is represented by a dot-dashed line 
in the left panel. A range of slit widths are displayed in each case: 1.0”, 0.67” & 0.33” (top to bottom respectively). 


signal to be detected in the simulated binary object, the pa¬ 
rameters used in these simulations are listed in Table 1. A 
percentage of the binary flux is allocated to a distortion in 
the PSF to simulate the effect of “uneven illumination”. 

Figure 6 shows an example result of the simulations 
when no distortion is present, the parameters for the posi¬ 
tion spectra shown are those given in Table 1, the specific 
values of Fdistort, ASdistor-t and seeing (FWHM) are 0.5, 
0.7” and 0.8” respectively. From this figure we can see that 
the “real” spectro-astrometric signal is unaffected by a re¬ 
duction in slit width, therefore the use of a narrow slit has a 
minimal effect on the extent of the real signal from the ob¬ 
ject and must therefore be an efficient method of reducing 
the artificial displacement in the observed position spectra. 

Also shown in Figure 6 is both the spectro-astrometric 
signal detected when the PSF suffered distortion, and the 
“real” spectro-astrometric signal that is detected when the 
PSF suffers no distortion (dot-dashed line). The simulations 
show that the extent of the detected artefacts is reduced by 
a factor of ~ 1.5 when the slit width is reduced by a third 
(from 1” to 0.67”), and when the slit width is reduced by 
two thirds (from 1” to 0.33”) the artefacts are reduced by 


a factor of ~ 3. Caveat: the use of a slit width smaller than 
the seeing will result in the loss of photons. As such, the 
measured positional accuracy is affected by the reduction in 
slit width (see Takami 2003). For example, if the slit width 
is reduced by a factor of 2 the number of photons detected 
(N) will reduce the signal to noise (oc W^) by a factor of 4, 
resulting in a reduction in the detection limit (oc 
and hence the positional accuracy achieved, by a factor of 
2 . 


4 CONCLUSIONS AND 
RECOMMENDATIONS 

Artificial spectro-astrometric signal can be produced in a 
position spectrum as a result of a distorted stellar image 
in a relatively wide slit, as was the case in data obtained 
using the 3.6m NTT. The observed artificial displacement 
looks similar to that produced when observing an object 
exhibiting bipolar structure, the only difference being that 
a real signal from bipolar structure will change sign when 
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observed using anti-parallel slit position angles whereas the 
artificial signal will not. 

Time variable distortion of the stellar image is a signif¬ 
icant problem, the time variation seen is presumably due to 
variation of the PSF, and as such the resulting magnitude of 
the effect is time-averaged. As a result the artefacts cannot 
be removed by subtracting position spectra obtained using 
long slit spectra with anti-parallel slit position angles, as is 
the case for other instrumental effects. Spectro-astrometric 
data must therefore be carefully monitored during observa¬ 
tions to determine when these artefacts arise and the obser¬ 
vations adjusted to compensate. 

Simulations show that the use of a slit width narrow in 
comparison to the seeing greatly reduces the extent of this 
instrumental effect. However, an image of high quality and 
good tracking of the telescope are mandatory in order to 
obtain position spectra with high positional accuracy. 

Based on observations and simulations we recommend 
that in order to minimise these artefacts and obtain position 
spectra with high positional accuracy the following are taken 
into account: 

(i) When observing use a slit width narrow in comparison 
to the seeing. This systematic effect was not observed in 
observations where the seeing is typically 1.5” - 2.0” and 
the slit width used was 1.0”. 

(ii) The use of an integral held unit (IFU) may be ad¬ 
vantageous as they do not experience uneven illumination 
(e.g. IFU’s with micro lenses or hbres), however although 
previous work by Garcia et al. (1999) appear to have been 
successful, instrumental effects arising in IFU’s may well be 
more complicated than those from long slit spectrographs. 
These should be investigated. 

(iii) Obtain spectra with a wide spectral range. The sys¬ 
tematic effect will create the same artihcial signal in a num¬ 
ber of emission and absorption features. Therefore if a wide 
spectral range is observed to include a number of absorp¬ 
tion and emission lines with a range of excitations it will be 
possible to deduce if the bipolar signal is real. 

(iv) Spectra should be obtained at anti-parallel slit po¬ 
sition angles. Not all instrumental effects are well under¬ 
stood or anticipated, as such comparison of position spectra 
obtained using anti-parallel slit position angles is the best 
method to determine whether or not the signal detected is 
due to real structure in the source. 
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APPENDIX A: PREVIOUS RESULTS 
REPORTED BY TAKAMI ET AL. (2001, 2003) 

The detection of bipolar outflows in Ha has been r e porte d 
previously in RU Lup and CS C ha. [ ^kami et alJ 11200111 : 
iTakami. Bailey & Chrysostomou! ll2003^ respectively, see 
Figure Al. In order to confirm that the results reported are 
due to real signals from the objects and not to the artefacts 
described in this work the (previously unpublished) anti¬ 
parallel slit data, which do not show the effect, are included 
in Figure Al. 

Displacement from the centroidal continuum position 
is seen in both objects at Ha A6563A, and additionally 
at [SH] A6716,6731A and [Nil] A6584A in RU Lup. Figure 
Al clearly shows that spectro-astrometric signal detected in 
anti-parallel slits (0° — 180°, 90° — 270°) are displaced in op¬ 
posing directions indicating a genuine signal from the object. 
In addition, the blueshifted Ha wing in RU Lup is displaced 
in the same direction as the forbidden emission, a well known 
probe for outflowing gas. The position angle of the displace¬ 
ment for Ha emission for CS Cha is perpendicular to the op¬ 
tical continuum polarisation, which is very often perpendicu- 
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Figure Al. Intensity and position spectra (Ho) of RU Lup (left) and CS Cha (right) reported in Takami (2001) and Takami (2003) 
respectively. Position spectra for each object were obtained at position angles of 0°, 90°, 180° &; 270°. 


lar to the jet axis iTakami. Bailev fc Chrvs ostomo^ ||200^). 
We thus conclude that the results reporte d bvlT akami et alJ 
^2001^ : lTakami. Bailev fc ChrvsostomotJ i200,'i[i ' are real and 
not artefacts. 


This paper has been typeset from a Tj^X/ f^Tj^X file prepared 
by the author. 
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ABSTRACT 

We demonstrate that artificial bipolar structure can be detected using spectro- 
astrometry when the point spread function (PSF) of a point source suffers distortion in 
a relatively wide slit. Spectro-Astrometry is a technique which allows us to probe the 
spatial structure of astronomical sources on milliarcseond (mas) scales making it pos¬ 
sible to detect close binaries and to study the geometry and kinematics of outflowing 
gas on scales much smaller than the seeing or the diffraction limit of the telescope. It 
is demonstrated that distortion of the PSF, caused by tracking errors of the telescope 
or unstable active optics during an exposure can induce artificial signals which may 
be misinterpreted as a real spectro-astrometric signal. Using simulations we show that 
these may be minimised by using a narrow slit relative to the seeing. Spectra should 
be obtained at anti-parallel slit position angles (e.g., 0° and 180°) for comparison in 
order to allow artificial signatures to be identified. 

Key words: Spectro-Astrometry Line: Formation Line: Profiles Stars: Pre-Main 
Sequence ISM: Jets and Outflows. 


1 INTRODUCTION 

Spatial resolutions of optical-IR facilities have dramatically 
improved over the last decade. Indeed, the Hubble Space 
Telescope and adaptive optics on ground based 8-lOm tele¬ 
scopes have provided spatial resolutions of 0.05 arcseconds, 
leading to a better understanding of the nature of a variety 
of astronomical objects. Despite this, even higher spatial 
resolutions are desired to study the geometry of extrasolar 
planetary systems, populations and formation of close bi¬ 
naries, mechanisms of mass ejection and accretion in young 
and evolved stars and the nature of active galactic nuclei. 
Optical-IR interferometry has begun to provide resolutions 
on milliarcsecond (mas) scales, however its applicability is 
limited due to the complexity of the technique. 

Spectro-astrometry is an alternative approach to study 
milliarcsecond structures at optical-IR wavelengths. The ba¬ 
sic concept is to measure the relative position of the source as 
a function of wavelength using either an intensity weighted 
centroiding algorithm or profile fitting (see Bailey 1998a for 
a more detailed description). The resultant “position spec¬ 
trum” will show structure on any spectral feature that is 
displaced from the centroid of the continuum source: e.g., 
an emission line arising in a binary companion, outflowing 
jets or any other structure that is not perfectly symmetric 
about the centroidal source, or whose displacement from the 

* E-mail: e.brannigan@star.herts.ac.uk 


source varies with wavelength. The smallest spatial scale ob¬ 
served is limited by the seeing or the diffraction limit of the 
telescope, however this technique indeed provides astronom¬ 
ically useful information at spatial scales much smaller than 
these. 

The concept was introduced in the 1980’s as “Differen¬ 
tial Speckle Interferometry” or “Chromatic Position Differ¬ 
ence” (e.g.,Beckers 1982), and demonstrated to detect close 
binaries (e.g., Sorokin & Tokovinin 1985). This broad band 
technique was limited by the necessity to correct for atmo¬ 
spheric dispersion which can be much larger than the ob¬ 
served structure. This difficulty is removed by measuring the 
displacement of a spectral line at high resolution. Aime et 
al. (1988) were able to measure spectro-astrometric struc¬ 
ture using special instrumentation. See Bailey (1998a) for 
review of these works. 

More recently Bailey (1998a) adapted this technique 
using a standard long-slit CCD spectrograph and coined 
the term “spectro-astrometry”. The method is similar to 
that used by Solf & Bohm (1993), who studied a jet from 
a young stellar object (YSO) on subarcsecond/arcsecond 
scales. Bailey (1998a) revised this method and has shown 
that with good spatial uniformity of the CCD coupled with 
a pixel scale that allows excellent sampling of the seing pro¬ 
file positional accuracies as small as Imas are achieved (see 
Takami 2003). Bailey has also shown that the technique is 
a powerful tool for discovering binary companions toward 
young and evolved stars, separating their spectra, and study- 
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ing the kinematics of YSO jets/winds and narrow line re¬ 
gions of active galactic nuclei (AGN). Since this method 
doesn’t require any special instrumentation, its popular¬ 
ity has steadily increased in recent years. Studies made for 
young binaries and YSO jets include Bailey(1998b), Takami 
et.al. (2001,2002,2003), Davis et al. (2001,2003), Whelan, 
Ray & Davis(2004) and Baines et al. (2004). Measurements 
using an integral field unit (IFU) were first conducted by 
Garcia, Thiebaut & Bacon (1999). A similar approach using 
VLTI-AMBER is investigated by Petrov (2003), Marconi et 
al. (2003) and Dougadous et al. (2003). Since the positional 
accuracy achieved depends on the PSF size as well as the 
photon noise, using an interferometer for this approach has 
the potential to achieve information on spatial structures at 
microarcsecond scales (Bailey 1998a). 

Instrumental effects which may compromise spectro- 
astrometric observations have not been fully investigated. 
These effects may indeed create false detections of close bi¬ 
naries or bipolar outflows. Possible instrumental effects in¬ 
clude; misalignment of the spectrum with CCD columns, any 
departure of the CCD pixels from a regular grid, imperfect 
fiatfielding or charge transfer deficiencies in the CCD (Bai¬ 
ley 1998b; Takami et al. 2001). To eliminate instrumental 
effects, Bailey (1998a,b) proposed to obtain position spec¬ 
tra at anti-parallel slit position angles (e.g.,0° and 180°) via 
rotation of the instrument. On comparison of these spectra, 
any real signal from the object will change sign, while the 
instrumental effects remain constant. We can thus eliminate 
instrumental effects by subtracting one position spectrum 
from the other. This method is followed by Takami et al. 
(2001,2002,2003) and Baines (2004). Garcia et al. (1999), 
Davis et al. (2001,2003) and Whelan et al. (2004) however 
eliminate instrumental effects by polynomial fitting of the 
continuum position. 

Bailey (1998a) reports a systematic effect upon obser¬ 
vation of sharp unresolved lines. This is caused by either 
telescope tracking errors or unstable active optics during 
the exposure, and appears in the position spectrum as a 
P Gygni type profile. Bailey warns that the target spectral 
lines should be well resolved to observe true positional dis¬ 
placement. This instrumental effect could be eliminated in 
some cases by Bailey’s method, but this cannot be done by 
polynomial fitting of the continuum position. In this paper 
we show that this instrumental effect can appear even if 
the line profile is fully resolved as was the case in position 
spectra of RU Lupi, obtained using the 3.5m New Technol¬ 
ogy Telescope (NTT). Takami et al. had previously detected 
the presence of bipolar outflow in this object using the 3.9m 
Anglo-Australian Telescope (AAT). Indeed the observed po¬ 
sitional displacement in the NTT looks similar to that seen 
in the AAT data. We show that the signal detected in the 
NTT data is not real, motivating us to investigate the effect 
of uneven illumination on spectro-astrometry, how to avoid 
it and whether previously detected bipolar outflows are in 
fact real. 

In §2 we describe the spectro-astrometric observations 
obtained using EMMI on the 3.6m NTT and the results 
in which artificial signatures were detected. These artefacts 
are identified by qualitatively inspecting position spectra ob¬ 
tained using anti-parallel slit position angles. In §3 we (1) 
show that the observed position spectra are explained by 
distortion/motion of the stellar image in a relatively wide 


slit, (2) show that the extent of the artefacts, in a given set 
of observing conditions, can be greatly reduced by using a 
slit width narrow in comparison to the seeing, and (3) show 
an example of the effect these artefacts have on real spectro- 
astrometric data. In §4 we give some recommendations on 
how to monitor the data for the false signatures of bipolar 
structure. In Appendix A, we show that the bipolar outflows 
reported by Takami et al. (2001,2003) are real and not due 
to the instrumental effect described here. 


2 OBSERVATIONS AND RESULTS 

Medium resolution spectroscopy of RU Lup was obtained 
with the 3.5m NTT in the early morning of April 8 
2004, when the seeing was measured to be 0.8”. The ob¬ 
servations were made using the REad Medium Disper¬ 
sion (REMD) mode of the ESO Multi-Mode Instrument 
(EMMI) with the ^6 grating and an arcsecond slit. The 
wavelength coverage spanned 6100 to 6800A with resolv¬ 
ing power = 5000. The spectra were obtained at 4 position 
angles: 0°, 90°, 180°, 270° by rotating the instrument. HeAr 
lamp spectra were used to calibrate the wavelength scale and 
a pixel scale of 0.165" with the MIT detector provided good 
sampling of the seeing profile. Flat fields were obtained by 
combining many exposures of the spectrograph illuminated 
by a halogen lamp. 

The data were reduced using the FIGARO package pro¬ 
vided by Starlink, following Bailey (1998a,b) and Takami 
(2001,2003). After subtracting the bias level and fiat- 
fielding, “position spectra” were determined by fitting the 
seeing profile at each wavelength with a gaussian function. 
Global instrumental effects in the position spectra are re¬ 
moved by polynomial fitting. This allows us to see the posi¬ 
tional displacement of any emission/absorption features at 
each slit angle relative to the continuum. In addition to the 
position spectra, the intensity spectra were obtained by sub¬ 
tracting the bias, fiat-fielding, subtracting the adjacent sky 
and extracting bright columns on the CCD. 

The spectrum of RU Lup shows a wealth of emis¬ 
sion lines as observed by Lago & Penston (1982) and 
Takami (2001). The intensity and position spectra of the 
Ha A6563A, Hel A6678A & [SII] A6731A emission lines and 
Lil A6708A absorption line, obtained using EMMI-NTT, are 
shown in Figure 1. The Ha emission line shows a large offset 
in the position spectrum, ~25mas from the continuum posi¬ 
tion and the blue and red shifted wings are displaced in dif¬ 
ferent directions. Similar displacement in Ha is reported in 
Takami (2001). Such structure usually indicates the presence 
of a bipolar outflow (see Appendix A). However a qualitative 
inspection of position spectra, obtained using anti-parallel 
slits (i.e. 0°&180°) in this data, reveals that the displace¬ 
ment seen in anti-parallel position spectra have equal sign 
indicating that the signature must be systematic. In addi¬ 
tion, the displacement observed with slit position angle 0° 
is less than that observed with position angle 180°. In fact 
two position spectra observed at a position angle of 0° ex¬ 
hibit dissimilar profiles indicating that the effect responsible 
for the spurious displacement signal is time variable and as 
such cannot be removed even by obtaining position spectra 
at anti-parallel slit angles. 

In addition to Ha A6563A and Hel A6678A, evidence for 
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Figure 1. Intensity and position spectra of RU Lup obtained using EMMI-NTT. From left to right; Top: Ha A6563A, Hel A6678A 
& [SII] A6731A emission and Lil A6708A absorption lines. The intensities shown are arbitrary; Lower four panels: spectro-astrometric 
signal (mas) for each corresponding spectral feature detected at 4 position angles 0°, 180°, 90°&270° in the spectrum of RU Lup. The 
bottom two panels show position spectra for anti-parallel slit angles 90° & 270° in each of the spectral features, similarly in the third 
and fourth panels from the bottom position spectra for anti-parallel slit angles 0° & 180° in each of the spectral features is shown. In the 
blue wing of [SII] A6731A emission feature the angular displacement is seen to change sign in corresponding anti-parallel slits, whereas 
displacement near zero velocity of the [SII] A6731A line (marked with a dot-dashed line), in the permitted emission lines and the Lil 


blueshifted forbidden line emission reverses sign at anti¬ 
parallel slit positions (0° vs. 180° ; 90° vs. 270°). This is best 
seen in the [SII] A6731A emission line in Figure 1. However 
a large, sharp displacement is observed close to zero veloc¬ 
ity in the position spectrum at both [SII] A6731A and [OI]^ 
A6300A emission lines. This displacement is observed in all 
the position spectra for [SII] A6731A (marked in Figure 1 
with a dot-dashed line) and [01] A6300A. This sharp dis¬ 
placement does not reverse sign in anti-parallel slit position 
angles and therefore this is likely to be due to instrumental 
effects and not to actual signal from the object. 


3 INDUCED SIGNAL DUE TO A DISTORTED 
STELLAR IMAGE IN A WIDE SLIT 

The instrumental effects described in §2 were not observed 
in previous observing runs at the AAT described in Takami 
et al. (2001,2003). There are a number of key differences be¬ 
tween these earlier observations and those reported in §2. 
Primarily, the typical seeing conditions differ dramatically 
(0.8" and 1.5"/2.0") while a wide slit of 1.0" was used in 

^ not shown as the effect is in this line is the same as that in 
the [SII] line, but the data at the wavelength of [OI] is of lower 
quality 


A6708A absorption line does not change sign in anti-parallel slits. 



Figure 2. Image capture in slit viewer 7/4/2004 using 3.5m NTT. 
Pixel scale: 0.4". Left and Centre: distorted image of the Classical 
T Tauri Star (CTTS) SZ68. Right: distorted image of the CTTS 
SZ69 


bipolar structure is seen in other bright permitted lines in¬ 
cluding Fell A6238/6248/6456/6516A. The Lil A6708A ab¬ 
sorption line exhibits displacement similar to that of Ha and 
the bright permitted emission lines, although in any given 
spectrum the direction of displacement in Li I absorption is 
of opposite sign to that of the emission lines. For example, 
at a position angle of 180° the blue and redshifted compo¬ 
nents of the bright permitted emission lines show positive 
and negative displacement respectively, whereas the Li I ab¬ 
sorption line shows negative and positive displacement in 
the blue and red wings respectively. 

Conversely, the direction of displacement due to 
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Figure 3. Top left: Intensity spectrum of RU Lup obtained on 8/4/2004 using EMMI on the 3.6m NTT. Lower left: simulated position 
spectra for a distorted PSF using 3 slit widths, 0.85", 0.5" & 0.17", flux ratio Fi : Fi = 0.1 and two point sources located off centre of the 
slit with identical spectra. Lower right: geometry used in the simulations to show that a distorted PSF can reproduce well the artefacts 
seen in the data. Top right: Close up view of the Hel, [SII] and Lil emission and absorption lines.Lower top right: positional displacement 
due to Hel, [SII] and Lil obtained from the simulations, where Fi : Fi = 0.1 and slit width = 0.5". The simulated displacement is similar 
to that seen in the data. 


all of the observations. These artefacts are seen in observa¬ 
tions when the seeing is smaller than the slit width used. 
Secondly, the stellar image suffered intermittent distortion 
during the later observations, as shown in Figure 2, due to 
an unstable active optics system. In §3.1 it is explained how 
these combined effects can induce an artificial signal in the 
position spectrum, and the extent of the artefacts that can 
be expected to occur for a given set of observing conditions is 
investigated. In §3.2 an example of the effect these artefacts 
can have on real spectro-astrometric signal is simulated. 

3.1 Concepts and Simulations 

When the intensity distribution from a given source is non- 
uniform across the slit the resultant spectrum can be slightly 
blue or redshifted. This blue/redshifted emission due to “un¬ 
even illumination” or “inaccurate centering”, (e.g., Bacciotti 
et al. 2002; Ardila et al. 2002) arises from an offset of the 
light source from the center of the slit, resulting in a devi¬ 
ation of the angle of incidence according to da = Ayjfcoi, 
where fcoi is the focal length of the collimator. According to 
the grating equation: 

mA = (sin a ± sin/3) (1) 


deviation in the angle of diffraction which in turn is trans¬ 
lated to a shift in the spectrum. Therefore, the shift in the 
spectrum, dA, is proportional to the position of the light 
source from the center of the slit axis Ay, dA = const, x 
Ay. This constant value can be readily determined from the 
spectral resolution of the instrument. For the REMD mode 
of EMMI-NTT a 1.0" slit provides R= 5000, corresponding 
to AA = I. 2 A giving us the relation: 

dA [A] = 1.2 X dj/ [arcsec] (2) 

As such, if a stellar image is distorted or elongated in the 
slit, perhaps due to tracking errors of the telescope or un¬ 
stable active optics, the spectrum close to the edges of the 
slit is slightly blue or redshifted. This geometry would be 
replicated by a bipolar outflow located within the slit. In 
order to determine the affect of these artefacts on spectro- 
astrometric observations simplified simulations of two point 
sources were carried out with their positions as depicted in 
Figure 3. Given that the point sources have spectra Fi/(A) 
and F 2 /(A), where Fi and F 2 are constant and /(A) is a 
normalised spectrum, the spectra are recorded at the detec¬ 
tor as Fi/(A + AA) and F 2 /(A — AA). Thus the centroidal 
position at each wavelength is given by: 


^ ^ ^ Fi/(A + 1.2Ay) - F 2 /(A - 1.2Ay) 
2 "" Fi/(A + 1.2Ay) + F 2 /(A - 1.2Ay) 


(3) 


where m is the diffracted order, A the diffracted wavelength, 
a is the angle of incidence wrt. the normal and /3 the an¬ 
gle of diffraction wrt. the normal, this leads to a subsequent 
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Parameter 

Description 

Distortion 

Binary 

Binary & Distortion 

PS 

Pixel Scale 

0.17 

0.17 

0.17 

^min 

start wavelength 

6100 

6100 

6100 

^max 

end wavelength 

6850 

6850 

6850 


Binary position angle 

- 

45° 

45° 

^^bin 

Binary separation 

- 

0.1" 

0.1" 

Fi 

Flux of primary binary component 

1.0 

3.0 

3.0 

F 2 

Flux of secondary binary component 

- 

1.0 

1.0 

"^distort 

Percentage of flux in distortion 

0.1 - 0.9 

- 

0.1 - 0.9 

^-^distort 

Distortion position angle 

135° 

- 

135° 

^^distort 

Distance between distortion and “real image” 

0.1" - 1.0" 

- 

0.1" - 1.0" 

FWHM 

Gaussian simulated seeing 

0.2" - 1.7" 

0.2" - 1.7" 

0.2" - 1.7" 

Aw 

Slit width 

0.1" - 1.4" 

0.1" - 1.4" 

0.1" - 1.4" 


Table 1. Parameters used in the simulations. 


To confirm that the artefacts seen in this data are a 
result of uneven illumination of the slit simulations of a dis¬ 
torted PSF were carried out and position spectra obtained 
using equation (3), and the geometry described in Figure 
3. The simulations reproduce a number of the artefacts ob¬ 
served in the data (Figure 1), specifically (see Figure 3): (1) 
evidence for bipolar structure exists in the bright permit¬ 
ted emission lines; (2) displacement seen at the Li I A6708A 
absorption line has opposite sign to that of the emission 
lines in a given spectrum; and (3) sharp displacement close 
to zero velocity is exhibited in the forbidden emission lines 
[Of] A6300A and [SII] A6731A. The extent of the displace¬ 
ment varies with the ratio Fi:F 2 , and in fact the relative 
displacement in the red and blue wings in the emission lines 
in a given position spectrum also depends upon this ratio. 
The extent of the displacement dramatically decreases as 
the slit width becomes smaller (see Figure 3). 

The point source geometry depicted in Figure 3 does 
not provide a good representation of observing conditions, 
as such the seeing conditions are simulated by convolving 
the two point sources by use of a gaussian function. The 
simulated intensity distribution is then described as: 

F{x,y) = Gi{x,y) + G 2 {x,y) (4) 


where: 

Gi{x,y) 

G2{x,y) 


2.773 

(FWHMYtt 

2.773 

(FWHMYtt 


„ - 2.773 

Fi e 


{FWHM)^ 


„ - 2.773 

F 2 e 


{FWHM)^ 


where FWHM is the full width half maximum of the seeing 
profile, (A xi, A yi) is the position of the primary point 
source, (A X 2 , A y 2 ) the position of the secondary point 
source, Fi and F 2 are the flux of primary and secondary 
point sources respectively. Combining equations (3) & (4) 
results in: 


//F(.,,)/(A + AA(y))d.d, 


The two point sources, having being allocated the rel¬ 
evant relative flux, are convolved by our simplified model 
of the seeing and subsequently shifted according to equa¬ 
tion 2. A position spectrum is then obtained as described 
by equation 5 

A simulated position spectrum (Figure 4), obtained as 
described above, shows displacement on a scale smaller than 
those artefacts obtained by the unconvolved point source. 



Wavelength (A) 

Figure 4. Position spectrum obtained as described in §3.1, where 
the parameters slit width, seeing ASdistort and Fdisiort are 1-0", 
0.8", 0.7" and 0.5 respectively. 


This is because much of the light falls close to and on the 
center of the slit and as such suffers less of a “shift”. In this 
instance a 1.0” slit is simulated with FWHM of 0.8”, as was 
the case for those observations shown in Figure 1. The two 
point sources were defined 0.7” apart, with equal flux (i.e. 
Fi:F 2 = 1) and both off center of the slit, as described in 
Figure 3. The angular scale of the displacement simulated 
for Ha emission (~ 25mas) compares well with the artefacts 
observed in the position spectra obtained at the NTT. 

In order to investigate the limitations imposed on 
spectro-astrometry due to systematic effects, the extent of 
the artefacts that can be expected when uneven illumination 
of the slit occurs (perhaps due to tracking errors or instru¬ 
mental concerns) for a given set of conditions is simulated. A 
single point source is allocated at the centre of the slit with 
flux Fi and subsequently convolved by use of a gaussian 
function. A distortion in the resulting FSF is simulated by 
“relocating” a percentage of the flux attributed to the point 
source (Fdistort) to a new position off centre of the slit. The 
results of these simulations are shown in Figure 6, in which 
the full extent of the artificial signal, measured from peak 
to trough, is shown for a range of seeing, slit width, distance 
between the “real image” and distortion (ASdistort) and the 
amount of flux attributed to the distortion (F^istoT-t), the val¬ 
ues of which are given in Table 1. The results are shown for 
an Fdistort of 0.1 and 1.0 for given slit widths and ASdistort- 
From this we can clearly see that for small ASdistort where 
the distortion in the FSF cannot be detected by eye artifi¬ 
cial spectro-astrometric signal is still generated on a scale 
~ 5mas in Ha, above the typical detection limit of Imas 
for previous spectro-astrometric observations (Takami 2001, 
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Figure 5. Simulated angular displacement in Ha, measured from peak to trough, for a range of seeing, slit width, distance between the 
“real image” and distortion (ASdjstort) ^nd the amount of flux attributed to the distortion which gives the ratio Fi:F 2 ), the 

values of which are given in Table 1. Top: Artificial signal generated for a given slit width (solid lines) for a range of seeing values, with 
Fi:F 2 = 1.0 (left) and Fi:F2=0.1 (right). Bottom: Artificial signal generated for a given separation between “real image” and distortion, 
ASdiatort (solid lines) for a range of seeing values, with Fi:F 2 = 1.0 (left) and Fi:F 2 = 0.1 (right). 


2003). In fact when the input parameters are similar to the 
conditions during earlier observations described in Takami 
2001, 2003 (i.e. seeing 1.5"coupled with a slit width of 1.0") 
the displacement in Ha shows an angular scale of ~7 - 8mas, 
similarly above the typical detection limit. Therefore it is im¬ 
perative to obtain data at anti-parallel slit angles in order to 
determine, by qualitative inspection, if any detected signal 
is in fact from a real structure in the object, even when the 
PSF appears symmetric. 

The model suggests that the effect can be minimised 
by defocusing the telescope. For example, Figure 5 shows 
that given a slit width of 1.0" an increase of 50% in the 
seeing parameter (from 0.8" to 1.2") results in a reduction 
of the angular scale of the simulated artefacts by a factor of 
2. However, the use of a slit width narrow in comparison to 
the seeing reduces the extent of the artefacts to a greater 
extent. From Figure 5 it can be seen that given a seeing of 


0.8" a decrease in the slit width by 50% (from 1.0" to 0.5") 
results in a reduction of the angular scale of the simulated 
artefacts by a factor of 3. As such it is clear that a reduction 
in the slit width is the best option in order to minimise the 
magnitude of these artefacts. 


3.2 Binary Simulations 

In order to determine (1) what, if any, effect the use of a 
relatively narrow slit might have on the detection of any 
real spectro-astrometric signal from the object, and (2) how 
this systematic effect will affect real signals from the object, 
the convolved model was adapted with the introduction of 
an artificial binary companion. The intensity spectrum al¬ 
located to the secondary binary component is the same as 
that in the primary with line to continuum ratio increased 
by a factor of 2 in order to allow a “real” spectro-astrometric 
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Figure 6. Binary simulations with (right) and without (left) a distorted PSF. The true displacement is represented by a dot-dashed line 
in the left panel. A range of slit widths are displayed in each case: 1.0", 0.67" & 0.33" (top to bottom respectively). 


signal to be detected in the simulated binary object, the pa¬ 
rameters used in these simulations are listed in Table 1. A 
percentage of the binary flux is allocated to a distortion in 
the PSF to simulate the effect of “uneven illumination”. 

Figure 6 shows an example result of the simulations 
when no distortion is present, the parameters for the posi¬ 
tion spectra shown are those given in Table 1, the specific 
values of Fdistort, ASdistort and seeing (FWHM) are 0.5, 
0.7” and 0.8” respectively. Prom this figure we can see that 
the “real” spectro-astrometric signal is unaffected by a re¬ 
duction in slit width, therefore the use of a narrow slit has a 
minimal effect on the extent of the real signal from the ob¬ 
ject and must therefore be an efficient method of reducing 
the artificial displacement in the observed position spectra. 

Also shown in Figure 6 is both the spectro-astrometric 
signal detected when the PSF suffered distortion, and the 
“real” spectro-astrometric signal that is detected when the 
PSF suffers no distortion (dot-dashed line). The simulations 
show that the extent of the detected artefacts is reduced by 
a factor of ~ 1.5 when the slit width is reduced by a third 
(from 1” to 0.67”), and when the slit width is reduced by 
two thirds (from 1” to 0.33”) the artefacts are reduced by 


a factor of ~ 3. Caveat: the use of a slit width smaller than 
the seeing will result in the loss of photons. As such, the 
measured positional accuracy is affected by the reduction in 
slit width (see Takami 2003). For example, if the slit width 
is reduced by a factor of 2 the number of photons detected 
(N) will reduce the signal to noise (cx: N'^) by a factor of 4, 
resulting in a reduction in the detection limit (cx: 
and hence the positional accuracy achieved, by a factor of 
2 . 


4 CONCLUSIONS AND 
RECOMMENDATIONS 

Artificial spectro-astrometric signal can be produced in a 
position spectrum as a result of a distorted stellar image 
in a relatively wide slit, as was the case in data obtained 
using the 3.6m NTT. The observed artificial displacement 
looks similar to that produced when observing an object 
exhibiting bipolar structure, the only difference being that 
a real signal from bipolar structure will change sign when 
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observed using anti-parallel slit position angles whereas the 
artificial signal will not. 

Time variable distortion of the stellar image is a signif¬ 
icant problem, the time variation seen is presumably due to 
variation of the PSF, and as such the resulting magnitude of 
the effect is time-averaged. As a result the artefacts cannot 
be removed by subtracting position spectra obtained using 
long slit spectra with anti-parallel slit position angles, as is 
the case for other instrumental effects. Spectro-astrometric 
data must therefore be carefully monitored during observa¬ 
tions to determine when these artefacts arise and the obser¬ 
vations adjusted to compensate. 

Simulations show that the use of a slit width narrow in 
comparison to the seeing greatly reduces the extent of this 
instrumental effect. However, an image of high quality and 
good tracking of the telescope are mandatory in order to 
obtain position spectra with high positional accuracy. 

Based on observations and simulations we recommend 
that in order to minimise these artefacts and obtain position 
spectra with high positional accuracy the following are taken 
into account: 

(i) When observing use a slit width narrow in comparison 
to the seeing. This systematic effect was not observed in 
observations where the seeing is typically 1.5" - 2.0" and 
the slit width used was 1.0". 

(ii) The use of an integral field unit (IFU) may be ad¬ 
vantageous as they do not experience uneven illumination 
(e.g. IFU’s with micro lenses or fibres), however although 
previous work by Garcia et al. (1999) appear to have been 
successful, instrumental effects arising in IFU’s may well be 
more complicated than those from long slit spectrographs. 
These should be investigated. 

(iii) Obtain spectra with a wide spectral range. The sys¬ 
tematic effect will create the same artificial signal in a num¬ 
ber of emission and absorption features. Therefore if a wide 
spectral range is observed to include a number of absorp¬ 
tion and emission lines with a range of excitations it will be 
possible to deduce if the bipolar signal is real. 

(iv) Spectra should be obtained at anti-parallel slit po¬ 
sition angles. Not all instrumental effects are well under¬ 
stood or anticipated, as such comparison of position spectra 
obtained using anti-parallel slit position angles is the best 
method to determine whether or not the signal detected is 
due to real structure in the source. 
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APPENDIX A: PREVIOUS RESULTS 
REPORTED BY TAKAMI ET AL. (2001, 2003) 

The detection of bipolar outflows in Ha has been reported 
previously in RU Lup and CS Cha, Takami et al. (2001); 
Takami, Bailey & Chrysostomou (2003) respectively, see 
Figure Al. In order to confirm that the results reported are 
due to real signals from the objects and not to the artefacts 
described in this work the (previously unpublished) anti¬ 
parallel slit data, which do not show the effect, are included 
in Figure Al. 

Displacement from the centroidal continuum position 
is seen in both objects at Ha A6563A, and additionally 
at [SH] A6716,6731A and [Nil] A6584A in RU Lup. Figure 
Al clearly shows that spectro-astrometric signal detected 
in anti-parallel slits (0° — 180°, 90° — 270°) are displaced 
in opposing directions indicating a genuine signal from the 
object. In addition, the blueshifted Ha wing in RU Lup is 
displaced in the same direction as the forbidden emission, a 
well known probe for outflowing gas. The position angle of 
the displacement for Ha emission for CS Cha is perpendicu¬ 
lar to the optical continuum polarisation, which is very often 
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Figure Al. Intensity and position spectra (Ha) of RU Lup (left) and CS Cha (right) reported in Takami (2001) and Takami (2003) 
respectively. Position spectra for each object were obtained at position angles of 0°, 90°, 180° & 270°. 


perpendicular to the jet axis (Takami, Bailey & Chrysosto- 
mou (2003)). We thus conclude that the results reported 
by Takami et al. (2001); Takami, Bailey & Chrysostomou 
(2003) are real and not artefacts. 

This paper has been typeset from a Tj^X/ file prepared 

by the author. 


© 2002 RAS, MNRAS 000, 1-9 





